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Curvaton model completed 
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In an inflationary cosmology, the observed primordial density perturbation could come from the 
quantum fluctuations of another light "curvaton" field, rather than the inflaton. In this case, it is 
essential that the curvaton decays, converting its perturbation to an adiabatic perturbation. For the 
first time, we consistently account for this decay in the simplest curvaton model V(a) — ^m^a^ and 
point out that it gives rise to an important logarithmic correction to the potential. Moreover, the 
potential will also receive a correction from the thermal bath. As a consequence, the dynamics of the 
curvaton are substantially altered compared to the original model in the majority of the parameter 
space. It will therefore be necessary to re-calculate all the predictions of the original model. 



We observe in the cosmic microwave background (and 
infer from large scale structure) an adiabatic, superhori- 
zon, almost scale-independent spectrum of primordial 
perturbations. A common hypothesis is that these per- 
turbations originated as quantum fluctuations of a light 
scalar field when the whole Universe was inside the hori- 
zon during inflation. This field could be the inflaton it- 
self, or another light field, called the curvaton a [H-Q, 
whose energy density pa- is subdominant during inflation. 
With the forthcoming Planck Surveyor Mission precision 
data it will be possible to test not only different models 
of inflation but also alternative models for the origin of 
the primordial perturbation. 

In the curvaton scenario inflation is driven by some 
field whose contribution to the curvature perturbation 
is negligible. There is also another field, the curvaton, 
which is very subdominant during inflation, with some 
initial field value cr*. During inflation the curvaton ac- 
quires perturbations from its quantum fluctuations. Af- 
ter inflation ends and the inflaton decays, the universe 
is filled with radiation. However, the curvaton perturba- 
tion continues evolving, with the curvaton field oscillat- 
ing in its potential until it too decays into radiation. At 
this time the curvaton perturbation is converted into an 
adiabatic curvature perturbation of the Universe. The 
amplitude of the final perturbation, which should match 
observations, depends on both how long the curvaton os- 
cillates before it decays, and on the shape of the potential. 
At the time of the decay, the curvaton may have become 
the dominant energy form or it may still be subdominant. 

The simplest curvaton model is given by the quadratic 
potential 



(1) 



where m is the curvaton mass. Such models have 
been extensively studied (see for example Q and ref- 
erences within). Models with renormalisable and non- 
renormalisable self-interacting terms in the Lagrangian 
have also been considered (see Q and references within 
for details). However, since the curvaton must decay, the 
curvaton must couple to other fields, and the couplings 



to other degrees of freedom will generate corrections to 
the potential, rendering the naive treatment of a purely 
quadratic potential inconsistent in large areas of the pa- 
rameter space. 

In this letter, we take the simplest realisation Eq. ([1]) 
of the curvaton model and consistently implement the 
requirement of curvaton decay. We consider the pertur- 
bative decay of the curvaton, although we note that the 
decay could instead occur via parametric decay 0, 
The curvaton, a scalar field, can be coupled in principle 
to either another scalar field </) or a fermion field V'- These 
interactions should be present in the bare Lagrangian of 
the curvaton, proportional to 

g^a^'f? or gaVii/' , 

with some coupling constants. The perturbative two 
body decay width to massless particles is then given by 
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where g is the coupling strength and 7 measures the num- 
ber of available decay channels, higher order effects and 
can account for many body decays (thus depends on the 
particle content). Unless there are a large number of de- 
grees of freedom present, 7 can be expected to be of 0(1), 
and thus we take 7 = 1 here for simplicity. 

The curvaton model now has two parameters, the mass 
TO and the coupling constant g (or F), and two initial con- 
ditions, and (T*. One of these parameters should be 
determined from the requirement that the curvaton pro- 
duces a curvature perturbation with the observed ampli- 
tude, C — 10~^. Furthermore as we will discuss below, 
the initial field value cr* can be constrained by probabilis- 
tic arguments. 

The existence of couplings to other fields implies quan- 
tum corrections of the well-known Coleman- Weinberg 
type for the effective potential of the curvaton^. The 



In principle, there will be gravitationally induced corrections to 
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FIG. 1: A schematic demonstration of the evolution of the 
curvaton energy density for the full potential Eq. JSjl (solid) and 
for the original model Eq. {TJ (dashed). In order to obtain the 
same value of C^, the full potential requires a much larger initial 
energy density than is obtained by considering only the original 
model. This difference can be several orders of magnitude. In the 
case of the full potential, the evolution begins as slow roll in Vt 
(oc a~^), followed by oscillations in Vt (oc a~^), oscillations in 
AV (oc a"**) and finally oscillations in Vq (oc a~^). 



effective one-loop potential for the model (with effective 
coupling g to effectively massless fields) is given by 



AV{a) = /3 
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where /? is some number depending on the number of 
fields (7 is coupled to, and fi is the renormalization point. 
A coupling to bosons gives a positive contribution to /3, 
while a coupling to fermions gives a negative contribu- 
tion. For now, we take /? = f for simplicity. 

After inflation ends, the background is dominated by 
radiation. Because of the thermal background there is 
also a thermal correction to the potential. To lowest 
order (for T ^ m) it is given by 



(4) 



where a again depends on the number of fields the curva- 
ton couples to; here we take a = 1. Right after infiation 
(assuming instant thermalization of the infiaton decay 
products) the temperature is given by T « (MpiJ,)^/^. 
Because of the thermal bath, the curvaton could also 



the potential. These will either be proportional to R (which is 
zero in a radiation dominated background, as considered here) 
or higher order operators suppressed by Mp. Therefore these 
corrections do not affects the dynamics after inflation. 



decay by absorption; however, there are kinematical con- 
straints that prohibit absorption if the curvaton couples 
to only one field (or, in case of fermions, if the ther- 
mal masses of left handed and right handed fermions are 
equal) [9|. Here we neglect absorption for simplicity. 
The total one-loop effective potential is thus 



AV{a) + VTicr,T) . 



(5) 



Hence the curvaton does not oscillate in a harmonic po- 
tential. This modifies the effective equation of state and 
the scaling law of the energy density of the curvaton, as is 
demonstrated in Fig. ([T]). The figure shows the evolution 
of Pa where the solid line is for the case where the curva- 
ton energy density is dominated in turn by the thermal 
correction Vt (in slow roll), Vt (oscillating), AV (oscil- 
lating) and finally oscillations in Vq and an instantaneous 
decay. In comparison, the evolution ignoring both decay 
and thermal corrections is shown with a dashed line. In 
order to match the observed curvature perturbation C, 
the curvaton energy density pa must have a certain value 
at the time of decay. Thus it is clear that the two sce- 
narios produce the same curvature perturbation ^ only 
with very different initial conditions. As a consequence 
the ramifications of the simplest curvaton model are not 
revealed by the study of the quadratic form Eq. ([1]), as 
we will now discuss. 

In order to further illustrate the importance of using 
the full effective potential Eq. ([5]), we consider the pa- 
rameter space in terms of the observables m and F. For 
each combination {m,T, H^^}, we choose the initial field 
value (T* such that the curvature perturbation matches 
the observed value C ~ 10^'"' where possible. It is given 
by 



C 



where 
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3p<T 



3pcr + 4:p^ 
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In our analytical estimates, we ignore the evolution 
of the factor (log [ga'^ / jj?) - 3/2) in Eq. (O since it 
changes only logarithmically, and replace with unity. (We 
have checked that its effect on the scaling law of the radi- 
ation is negligible.) We also assume sudden decay when 
H = V and that only one of the terms in the potential 
dominates at a time. The curvaton is assumed to be 
frozen ct = cr» until = V" , where denotes the end 
of inflation throughout this letter. We further approxi- 
mate Tdec as Pa/4M^H^. 

To ensure a consistent scenario, we apply various con- 
straints to the parameter space: 

(i) the initial field value should be sub-Plankian (cr* < 



3 



10-' 



> 



10-^1 



WIMP 




10' 



10-- 



> 

o 



10-- 



10-' 



10-- 



10-- 



10' 



10' 



10- 



F-P 




10" 



10* 
/tttr / GeV 

(a) This figure is plotted with H, = GeV. The WIMP 
limit (condition (v) in the text) requires F to be above the 
dashed line. 
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(b) This figure is plotted with H, = 10^^ GeV. Note the 
different scale of the axes compared to (a). 



FIG. 2: The coloured region shows the allowed parameter space. If the contributions of AV and Vr were neglected, the allowed region 
would extend to include the area bordered by the dotted line. In (b) this extends beyond the scope of the figure. Below the solid line 
(marked "F-P") is the region with high probability (assuming a long period of inflation). 



(ii) the curvaton should be hght during inflation 
(TOe// < i?*), 

(iii) the curvaton energy density should be initially neg- 
ligible (r* < 1), 

(iv) the curvaton should decay before nucleosynthesis 
(r > 10^22 Qgv) 

(v) the curvaton should decay before dark matter 
freezeout. If dark matter is a WIMP which decou- 
ples at Tdm ~ 0(10) GeV, then T > IQ-^^ GeV. 
This bound depends on the model of dark matter 
and could be relaxed to allow smaller F. 

(vi) non-Gaussianity should be below the current obser- 



vational bound, usually parameterised by / 



NL 



5/(4rdec) < 74 [10|. 



(vii) the decay width should be at least as large as the 
gravitational decay width (F > \m? /Ad^ where we 
conservatively take A = lO^'^) 

(viii) the model should generate a curvature perturbation 
which matches observations {(, ^ 10^''^) 

Constraints (i) and (iii) do not impact the parameter 
space as the other constraints are stronger. 



If inflation lasts long enough, we can also constrain the 
initial condition for the curvaton from the Fokker-Planck 
limit. It can be shown that in the limit of an infinite 
duration of inflation, the initial field value will be limited 
by the equilibrium probability distribution 

8^2 
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which is the equilibrium solution of a Fokker-Planck 
equation balancing the classical dynamics and quantum 
noise. This is valid under the assumption that the Hubble 
parameter evolves slowly in the last 60 or so e-foldings 
of inflation. We consider the region P > to 
have high probability under the above assumption. To 
be precise, for a finite number of c-folds, one should also 
consider the rate of relaxation into the asymptotic dis- 
tribution Eq. ([5]); we do not address this issue here^. 

In Fig. ([2]), we show the available parameter space for 
two representative values of i?*: = 10^° GeV and 



^ The potential in Eq. ^ does in principle contain gravitationally 
induced quantum corrections. For simplicity, we do not consider 
these here and the Fokker Planck limit should be taken as in- 
dicative. 
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i?* = 10^^ GeV. The results are qualitatively the same 
even if values of a and /3 as large as 100 are used (the 
figures show a — {3 — 1). The coloured regions show the 
allowed parameter space for the model. This is bounded 
from above by the requirement of generating the observed 
C (viii) and from below by the observational bound on 
Inl (vi). In Fig. ((2a|) . the bound for large m is the 
gravitational decay width (vii) and in Fig. (|2b|) it is the 
lightness of the curvaton (ii). The low m bound is the 
most interesting and is due to Vt and AV. Below this 
lower bound, these terms dominate the potential initially. 
However, in this region it is not possible to generate a 
sufficient amplitude of perturbations while satisfying all 
other constraints. If the terms Vt and AV were ignored, 
the allowed region would fill the area outlined by the dot- 
ted lines. Thus, using the simplest, complete curvaton 
potential gives a substantially different parameter space 
than the assumption of a pure quadratic potential. 

We also show the region of the parameter space which 
has high probability (assuming a long enough period of 
inflation). This is below the solid line marked "F-P". 
This high probability region is greatly reduced from the 
parameter space allowed by the constraints (i) to (viii) 
above and favours low m and low F. 

The dynamics of the curvaton in the allowed region can 
also be different from the original model. Some of the low 
m allowed region is initially dominated by Vt- In addi- 
tion, AV may also play a role in the dynamics, although 
it is subdominant. This would affect the calculation of 
the non-Gaussianity parameters /^vl and gjvL- In a later 
paper, we expect to produce a full numerical analysis of 
the parameter space, including predictions for /atl and 
gNL- We expect to be able to carefully explore the whole 
curvaton parameter space, including all allowed values of 
i?* • We will also consider the corrections to other curva- 
ton models. 

In summary, we have completed the simplest curvaton 
model by including the full effect of the required decay 
of the curvaton. This gives rise to an additional term 
in the effective potential of the curvaton, proportional to 

log a. A thermal correction to the potential should 
also be considered. Using an analytical approximation, 
we were able to calculate the parameter space for this 
model consistent with various constraints, such as the 
observational limit on /jvl ■ We showed that the parame- 
ter space calculated for the original curvaton model Vq is 
greatly reduced when the thermal correction and the cor- 
rection AV due to the curvaton's coupling are included. 



Computations of curvaton parameters using only Vq are 
not necessarily valid due to the differing dynamics of the 
curvaton model. They should therefore be repeated with 
the realistic potential presented in this letter. 

A particularly large effect on the available parameter 
space comes from consideration of the initial conditions 
for the curvaton's evolution. Assuming a long enough 
period of inflation, we calculated the probability of the 
curvaton having the correct initial conditions in order to 
match the observed curvature perturbation. The high 
probability region is substantially reduced compared to 
the parameter space neglecting to consider the initial con- 
ditions. In conclusion, it is essential to consider the full 
potential for the curvaton as this can have a dramatic 
effect on the predictions of the model. 

OT acknowledges the support by the Magnus Ehrn- 
rooth foundation. RL and KE are respectively supported 
by the Academy of Finland grants 131454 and 218322. 
We wish to thank Tomo Takahashi for useful comments 
on an earlier version of this manuscript. 



t 
t 

[1 
[2: 

[s: 

[4: 

[5: 
[6: 

[7 

[8 
[9 

[10: 
[11 

[12 



Electronic address 
Electronic address 
Electronic address 



ihelsinki.fil 
ihelsinki.fi' 



(2002) 395- 
(2002) 5-14. 
(2001) 215- 



kari.enqvist(! 
rose.Iernerif 
olli .taanila@helsinki.fi 
K. Enqvist, M. S. Sloth, Nucl. Phys. B626 
409. hep-ph/0109214. 

D. H. Lyth, D. Wands, Phys. Lett. B524 
hep-ph/0110002 . 

T. Moroi, T. Takahashi, Phys. Lett. B522 
221. hep-ph/0110096 . 

K. Dimopoulos, G. Lazarides, D. Lyth, R. Ruiz de Austri, 

Phys. Rev. D68 (2003) 123515. hep-ph/0308015|. 

K. Enqvist, arXiv:1012.1711 [astro-ph.CO]]. 

K. Enqvist, S. Nurmi, G. I. Rigopoulos, JCAP 0810 

(2008) 013. arXiv:0807.0382 [astro-ph]]. 

M. Bastero-Gil, V. Di Clemente, S. F. King, Phys. Rev. 

D70 (2004) 023501. hep-ph/0311237 . 

A. D. Linde, Phys. Lett. B114 (1982) 431. 

See e.g. H. A. Weldon, Phys.Rev.D26:1394,1982; 

Phys.Rev.D26:2789,1982; P. Elmfors, K. Enqvist and L 

Vilja, Nucl. Phys. B421 (1994) 459. 

E. Komatsu et al. [ WMAP Collaboration ], Astrophys. J. 
Suppl. 192 (2011) 18. [arXiv:1001.4"538i [astro-ph.CO]]. 
K. Enqvist, S. Nurmi, G. Rigopoulos, O. Taanila, 
T. Takahashi, JCAP 0911 (2009) 003. arXiv:0906.3126 
[astro-ph.CO]]. 

A. A. Starobinsky, J. Yokoyama, Phys. Rev. D50 (1994) 
6357-6368. [astro-ph79407016]. 



